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ABSTRACT 

The detection of exoplanets through direct imaging has produced numerous new positive identifica- 
tions in recent years. The technique is biased towards planets at wide separations due to the difiicuhy 
in removing the steUar signature at smah angular separations. Planets in eccentric orbits will thus 
move in and out of the detectable region around a star as a function of time. Here we use the known 
diversity of orbital eccentricities to determine the range of orbits which may lie beneath the detection 
threshold of current surveys. We quantify the percentage of the orbit which yields a detectable sig- 
nature as a function of semi-major axis, eccentricity, and orbital inclination and estimate the fraction 
of planets which likely remain hidden by the fiux of the host star. 

Subject headings: planetary systems - techniques: high angular resolution - radial velocities 
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1. INTRODUCTION 

The diversity of exoplanetary systems has been re- 
vealed through a variety of detection techniques which 
are sensitive to different kinds of planets. The majority of 
planets have been detected indirectly via the radial veloc- 
ity (RV) and transit methods, but the technique of direct 
imaging has made important contributions by detecting 
young planets at wide separations (> 10 AU). Direct 
imaging is a difficult technique to execute and requires 
the use of adaptive optics and high contrast imaging to 
disti nguish the flux of the planet fr om that of the host 
star (|Oppenheimer fc Hinklevll2009| ). The characteristic 
properties of atmospherically-induced speckl e noise can 
be treated using adaptive opti cs techniques (.Racine et al.l 
119991 iMacintosh et al.l I2007D . Speckle noise due to 
telescope/optics imperfections require more advanced 
techniques such as angu l ar differential imaging (ADI) 
(jMarois et al.ll2005L [20061 : iLafren iere et al."2007a'). spec- 
tral differential im aging (SDI) (Lafrcnicre et al. 2007b; 
iVigan et al. I2010f). and polar imetric differential imag- 



ing (POn (jOuanz et al1l2012[ ). The list of both tech- 
nical and scientific achievements from exoplanet imag- 
ing surveys is far too exhaustive to describe here in 
detail, but include such milestone exoplanet discover- 
ies a s beta Pic b ([Lagrange et al.l l2010t iCurrie et al] 
201 ip the HR 8799 planetarv svstem (iMarois et al.ll2008r 
20ld iHinklev et al.l 120111: ICurrie et al.ll2012aD. a nd Fo- 



malhaut b ( Kalas et al.l 120081 : ICurrie et al 

The star-planet separations explored by the imaging 
and RV techniques are complementary to each other 
since the RV method is biased towards shorter orbital 
periods. Although the imaging technique requires con- 
siderable follow-up observations to constrain the or- 
bits of the planets detected, RV planets provide com- 
plete Keplerian orbital solutions at the time of discov- 
ery which may be used to study the overall statistics 
of exoplanetary systems. The eccentricity distribution 
has b een investig ated for RV planets (,Shen fc Turneii 
120081 [Hogg et alj 120101) and also compared to that of 
the K epler candidates ([Moorhead et aT1l2011l : [Kane et ahl 
[20T1 . The origin of the observed distribution has 
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also been studied in the context of plan et forma- 
tion scenarios and early dynamical evolution ([Namounil 
12005"; 'Ribas & Miralda-Escude"2007: Ford & Rasio|[2008t 
Malnibcrg & Davics 2009; Kley & Nelson 201^n:f this 
eccentricity distribution applies also to the kinds of plan- 
ets detected via direct imaging then the detectability of 
those planets will depend on where they are in their or- 
bit since the star-planet separation is time dependent 
(Brandckcr et al. 2006; Bonavita et al. 20i3)- 

Here we investigate the effects of exoplanet orbital ec- 
centricity on the detection efficiency of imaging surveys. 
We describe the constraints on observing planets in ec- 
centric orbits depending upon the sensitivity of direct 
imaging experiments to the inner region of the system 
and the cadence of the observations. We calculate the 
percentage of the orbit for which the planet will be de- 
tectable as a function of eccentricity and the sensitivity 
threshold. We study the known eccentricity distribution 
as a function of semi-major axis and star-planet angular 
separation. Finally, we use this distribution to estimate 
the detection efficiency of direct imaging for eccentric 
exoplanets. 

2. OBSERVATIONAL CONSTRAINTS 

In addition to the science requirements of imaging ex- 
periments for achieving high-contrast observations that 
allow one to detect the extreme star-planet flux ratio, 
one also needs to contend with the limit of how close to 
the host star can reasonably be probed for planets. The 
angular resolution needed for achieving this is limited 
by the diffraction-limited resolution where the objects 
are separated by X/D where D is the aperture diame- 
ter of the telescope. In practice, there are further tech- 
niques that may be brought to bear on the problem that 
can further suppress the starlight and enhance the con- 
trast ratio. These techniques vary enormouslessly in both 
their effectiveness and their number. Here we quantify 
the problem in terms of a uniform inner region around 
a star for which planet detections are inaccessible. We 
refer to this region as an exclusion zone with radius r^- 
Although a practical exclusion zone will be quite compli- 
cated non-uniform pattern due to diffraction and speckle 
effects, this allows us specify a lower-limit which may be 
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Fig. 1. — Three example orbits for a planet in orbit around a 
direct imaging target. The shaded region indicates the angular re- 
gion around the star for which the example imaging experiment 
is insensitive to positive planet detections (exclusion zone). The 
orbits each have a semi-major axis equal to the radius of the exclu- 
sion zone and have eccentricites of 0.3, 0.6, and 0.9. The orbits are 
shown as points which indicate the planet location in equal time 
increments. The numbers in parentheses show the perecentage of 
the orbit that lies outside the exclusion zone. 

adapted to any particular experiment. 

The problem of detecting eccentric exoplanets for 
which the star-planet separation is time-dependent is 
portrayed in Figure [1] The shaded region represents 
the exclusion zone for a given star. The three orbits 
shown represent three Keplerian orbits which have identi- 
cal components with the exception of eccentricity, shown 
here for e = 0.3, 0.6, and 0.9. The argument of perias- 
tron is set to a; = 315° for this example. The semi-major 
axis in all cases is set to a = re, such that a circular 
orbit would lie on the threshold of detectability. The 
orbits are shown as points which indicate the planet lo- 
cation in equal time increments to illustrate the relative 
amount of time the planet spends at different phases of 
its eccentric orbit. In this example, there are clear ad- 
vantages for highly eccentric planets whose orbits result 
in a larger angular separation from the host star and a 
larger percentage of the orbit outside of the exclusion 
zone. 

The percentage of the orbit during which the planet 
is in the detectable region depends on the eccentricity 
and semi- major axis relative to r^. This is shown in 
Figure [2] for four different values of a (in units of Tg). For 
a < there are large swathes of eccentricities for which 
the orbit never exits the exclusion zone. The limit at 
which the planet never exits the exclusion zone is found 
from evaluating the star-planet separation at apastron, 
that is where a(l -f e) < rg. As e approaches unity, 
this limit is located at a = 0.5re. A similar bound on a 
semi-major axis beyond which the planet never enters the 
exclusion zone is more open-ended since the periastron 
boundary condition of a(l — e) > approaches infinity 
as e approaches unity. 

Another observation constraint is that of the observ- 




Eccentricity 

Fig. 2. — The dependence of the percentage of an orbit spent 
outside the exclusion zone on the orbital eccentricity for four dif- 
ferent semi-major axis. The semi-major axes are in units of the 
exclusion zone radius, r^. 

ing cadence. The planet is moving at its slowest speed 
when the star-planet separation is approaching maxi- 
mum. Thus even sparse sampling should in principle 
allow the tracking of the planetary orbit. There are how- 
ever two things to consider here. Firstly, if the movement 
of the planet is less than the positional accuracy of the 
planetary point-spread-function then this will create de- 
generacy in the derived Keplerian solution for the planet. 
On the other hand, too sparse sampling could result in 
the planet being missed in future (or past) observations 
if the apastron is near the detection threshold, such as 
the a — O.Gr-e case shown in Figure [2l Note that the ex- 
amples shown here aid in defining the problem posed by 
eccentric orbits, but assume face-on orbits for the specific 
numbers calculated. We discuss the effects of inclination 
in Sectional 

3. ECCENTRICITY DISTRIBUTION AT LONG PERIODS 

Here we briefly discuss the eccentricity distribution of 
the known exoplanets, particularly at relatively long pe- 
riods. To do this, we utilize the orbi tal parameters store d 
in the Exoplanet Data ExploreiQ (| Wright et all 120111) . 
We extracted the data of 513 planets, along with the host 
star properties, with data current as of 20th November 
2012. These data include all planets with the necessary 
Keplerian orbital solutions as well as stellar distances. 

In Figure [3] we have plotted these data to show the dis- 
tributions of eccentricity values. The left plot shows this 
distribution as a function of semi-major axis, where we 
have logarithmically scaled the size of the points depend- 
ing on the planetary mass. The relative dearth of planets 
between 0.1 AU and 1 .0 AU is well-known, noted for ex- 
ample bv lButler et al.l (|2006,). As is clear from the plot, 
the mean eccentricity increases with semi-major axis and 
is 0.06, 0.24, and 0.27 for semi-major axis ranges of 0.0- 
0.1, 0.1-1.0, and 1.0-2.0 respectively. What is not so 
clear is the increase of planet mass with eccentricity; 2.5, 
3.6, and 5.4 Jupiter masses for eccentricity ranges of 0.0- 
0.3, 0.3-0.6, and 0.6-1.0 respectively. This is of particu- 
lar relevance to the imaging surveys since they are more 
sensitive to more massive planetary companions. 

^ |http : //exoplanets ■ org/] 
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Fig. 3. — Left panel: The eccentricity of the known exoplanets as a function of semi-major axis. Right panel: The angular separation 
(arcsecs) of the same sample of planets at apastron and assuming a face-on orbit (i = 0°). In each plot, the size of each plotted point is 
logarithmically proportional to the mass of the planet. 



In the right panel of Figure [3] we have plotted the ec- 
centricity against the calculated angular projected sepa- 
ration of the planet from the star. This angular separa- 
tion is calculated for a face-on orbit when the planet is 
located at maximum separation (apastron). For a Kep- 
lerian orbit, the time-dependent star-planet separation, 
r, is given by 



where / is the true anomaly. The general expression (in 
radians) for the angular separation as a function of time 
is then 

AO =-( cos^ (w + /) + sin^ + f) cos^ i)^ (2) 

where uj is the periastron argument, i is the or- 
bital inclination, and d is the star-observer distance 
(|Kane fc Gelinoir201lD . This angular separation is a use- 
ful diagnostic to show which of the known exoplanets 
may yield positive results from a certain imaging experi- 
ment. What is clear is that many of those planets which 
fall into a detectable range of angular separations will be 
at the apastron of an eccentric orbit. 

4. DETECTION EFFICIENCY FOR ECCENTRIC ORBITS 

We now utilize the components of the previous sections 
to calculate the expected fraction of planets that will be 
detectable from images experiments. To do this, we re- 
quire the inclusion of the effects of inclination since this 
can have a large effect on the projected separation be- 
tween the star and planet. Calculations from Equation [2] 
are combined with those of a Keplerian orbital solution 
to compute the projected separation as a function of or- 
bital phase, where phase zero is the location in the orbit 
where lu + f — 270°. We also set w = 315° and a = r^. 
The choice w is to demonstrate the asymmetry that can 
result when considering eccentric orbits, a more thorough 
discussion with application to photometric phase curves 
of which may be found in iKane fc Geli^ (l20Tll ). We 
show the results of these calculations in Figure |3] where 
we use three different values for eccentricity (e = 0.3, 
0.6, and 0.9) and inclination {i = 0° (face-on), 45°, and 



90° (edge-on)). The shaded region indicates the exclu- 
sion zone (as for Figure [1]). The eccentricity, inclination, 
and percentage of the orbit that lies outside the exclusion 
zone, respectively, are given at the top of each subplot. 
In general, the most favorable condition for detection is 
that for which inclinations are close to face-on, although 
high-eccentricity orbits suffer less in this regard when 
orbits become close to edge-on. Depending on the asym- 
metry in the projected separation (sensitively dependent 
on w), increased values of a will result in multiple oppor- 
tunities for detection as the orbit moves in and out of 
the exclusion zone. 

To quantify the effect of all these aspects on the de- 
tection efficiency of imaging experiments, we performed 
a detail Monte-Carlo simulation. This simulation calcu- 
lates the percentage of the orbit which lies outside the 
exclusion zone for each of the known exoplanets described 
in Section [3l We calculate this a function of the exclu- 
sion zone expressed as an angular separation which uses 
the known distance to the star. We also consider three 
inclinations for all planets; i — 0°, 45°, and 90°. We 
then randomly determine if the planet lies outside the 
exclusion zone, thus using the percentage of time outside 
the exclusion zone as a probability distribution. The re- 
sults of this simulation are shown in Figure [5] where we 
have calculated the percentage of the planets recovered 
in each of the simulations for each exclusion radius. We 
consider the range of exclusion radii from 0.01"to 0.2", 
a range for which many of the known exoplanets, par- 
ticularly those at short periods, will remain elusive. For 
comparison, the angul ar separation of H R 8799 e from 
the host star is 0.37" (|"Marois et al.ll2010( ). and the inner 
working angle of selected future ground-ba.sed i maging 
instruments described by iBeichman et al.l ()2010D range 
from 0.03"to 0.17". The results of this simulation sup- 
port the conclusion from Figure |4] that face-on inclina- 
tions are more favorable for detection. More importantly 
though, it demonstrates the dramatic dependence on the 
exclusion radius for the known exoplanets. This depen- 
dence on the exclusion radius is true even if all the orbits 
are circular. Note that assuming an isotropic distribu- 
tion of orbital inclinations by randomly selecting the in- 
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Fig. 4. — The projected star— planet separation as a function of orbital phase for eccentricities of e = 0.3, 0.6, and 0.9 and inclinations of 
i = 0° (face-on), 45°, and 90° (edge-on). We use ui = 315° and a = re in this example. The shaded region indicates the exclusion zone (as 
for Figure [T] and the right-hand number in each plot is the percentage of the orbit that lies outside the exclusion zone. 

We first used the eccentricities of the planets between 
0.0-0.1 AU to represent near circular orbits. This con- 
sists of 165 planets which we use as the seed for the 
Monte-Carlo simulation whilst fixing all semi-major axes 
to a given value. As for the previous simulation, we 
choose random points in the orbital phase for each planet 
and determine if the planet is detected (outside the ex- 
clusion radius) or undetected (inside the exclusion ra- 
dius). For semi- major axes of a = l.Org, 76% of planets 
are detected, whereas for a = 1.5re, 100% of planets 
are detected. We then repeated this simulation by using 
the eccentricities of the planets which lie between 1.0- 
2.0 AU; a sample which is more representative of planets 
at longer orbital periods where imaging experiments are 
more sensitive. For semi- major axes of a = l.Org, 62% 
of these planets are detected, and for a = 1.5re, 96% of 
planets are detected. This is consistent with the results 
shown in Figure [2] and indicates that the vast majority 
of eccentric planets with a > l.Srg should be detectable 
at any given time. Note that this assumes that the ec- 
centricity distribution between 1.0-2.0 AU continues at 
longer orbital periods. Figure [3] suggests that the mean 
eccentricity may decline beyond 2.0 AU. However, it is 
unclear if this is an astrophysical reality or a symptom 
of RV su rvey incompleteness to eccentric orbits at longer 
periods ()Cumminell200l IShen fc Turneri[2008i ). 




Exclusion Radius (") 

Fig. 5. — The results of a Monte-Carlo simulation which calcu- 
lates the percentage of the known exoplanets (see Section [3]l which 
would be detected using direct imaging as a function of the exclu- 
sion radius for that particular experiment. Since we cannot assume 
an inclination, we have performed this simulation for inclinations 
of i = 0°, 45°, and 90°. 

clination for each planet in the Monte-Carlo simulation 
results in a dependence that is roughly equivalent to the 
i = 45° case shown in Figure [S] 

Finally, we performed a simulation which uses the ob- 
served eccentricity distribution described in Section [3l 
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5. DISCUSSION 



There are numerous imaging surveys which have con- 
sidered the detectio n Umitations of their particular ex- 
perim e ntal designs (iChauyin et al.ll2010l: iBonavita et alj 
2012t iBrandeker et alJ l2006t iDelorme et all 120121: 



Janson et al. Il 2012t Lafreniere et al.|[2007d iNielsen et alJ 
2008; Nielsen fc Closd 120101 : IVigan et all I2012n . Th"ese 

often treat the eccentricity distribution as a basic 
gaussian function based upon the RV exoplanets. 
Examples of the methods explored are the commonly 
used probability dens ity distribution a, p proac h described 
in Appendix A of IBrandeker et al.l (|2006) and the 
Monte Carlo simula tion code (MESS) developed by 
IBonavita etHI ()2012( ). The differ ence with our approach 
here is that (a) we use the RV eccentricity distribution 
directly rather than infer a function from that distri- 
bution, and (b) we specifically determine the effect of 
eccentricity and the correlated parameters of periastron 
argument and inclination on detection efficiency. 

The Monte-Carlo simulation we perform in Section [4] 
uses the eccentricity distribution of RV planets at long- 
periods as a proxy for the expected distribution at even 
larger separations. However, as was noted at the conclu- 
sion of the Section, it is unclear if this distribution does 
indeed persist in this fashion. Additionally, the sam- 
ples of RV host stars and imaging survey targets tend to 
represent quite different stellar demographics in terms 
of their dista nces and ages. Ecceii t ricity studies of bi- 
nary stars bv lDuauennov fc Mavoii ()199lD indicate that 
the distribution peaks at e ~ 0.35 for P > 1000 days, 
whereas the mean eccentricity of the known RV planets 
between 1.0-2.0 AU is 0.27 (see Section jS]). However, 
considering the very different formation mechanisms be- 
tween stars and planets, such a comparison may be of 
dubious value. Cons traints on the o rbital eccentricity 
of beta Pictoris b by iChauvin et al.l ()2012D and on the 
eccentricities of the HR 8799 planets by iCurrie et al.l 
(j2012a) have shown that these eccentricities are rela- 
tively low (< 0.3). Thus the number of imaged plan- 
ets are small enough and the eccentricities low enough 
such that little may be discerned regarding the eccentric- 
ity distribution from these planets alone. Increasing this 
sample by a factor of 10 would start to produce meaning- 
ful insight into the distribution of eccentricities beyond 
that which is currently being probed by RV surveys. As 
imaging experiments are able to decrease the exclusion 
radius and provide adequate observational cadence, the 
ability to detect the orbital motion of detected planets 



(and thus constrain the Keplerian properties of that or- 
bit) will correspondingly increase. 

6. CONCLUSIONS 

The exoplanet detection method of direct imaging is a 
powerful technique which is rapidly developing in both 
technical and scientific contexts. The major frontier push 
is to not only to detect smaller masses, but to detect 
those at smaller separations from their stars. Thus the 
overlap between the semi-major axis parameter-space of 
the RV and imaging techniques is becoming larger at a 
rate such that it's prudent to investigate the kinds of 
Keplerian orbital solutions the imaging technique may 
expect to encounter with new discoveries. Here we have 
shown the impact of orbital eccentricity on exoplanet de- 
tectability through simulations, the observed distribution 
of Keplerian orbital parameters, and the use of a simplis- 
tic model of the exclusion zone for a particular experi- 
ment. In practice, the exclusion zone is not well-defined 
and is a complex function of the diffraction and speckle 
effects in the resulting images. The exclusion zone used 
here though well represents an inner region where detec- 
tions are rendered virtually impossible by the flux of the 
star. These results show that the eccentricity can have 
a profound effect of the detection of eccentric exoplanets 
whose semi-major axis lie in the range 0.5 < < 1.5. 
We have also quantified the effects of inclination on these 
results and shown how face-on orbits can increase the de- 
tection efficiency for eccentric orbits. As the techniques 
for direct imaging improve and further epochs are ac- 
quired for the known imaged systems, the precise effect 
of eccentricity on these experiments will be become abun- 
dantly clearer and add enormously to our knowledge of 
the eccentricity distribution as very long orbital periods. 
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